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SUMMARY 


The governing equations for an unsteady turbulent boundary layer on a swept infinite 
cylinder, composed of a continuity equation, a pair of momentum equations and a pair 
of turbulent energy equations which include upstream history efforts, are solved 
numerically. An explicit finite difference analog to the partial differential equations 
is formulated and developed into a computer program. Calculations were made for a 
variety of unsteady flows in both two and three .li;. - r 'sions but primarily for two 
dimensional flow fields in order to first understand some of the fundamental physical 
aspects of i-.isteady turbulent boundary layers. Oscillating free stream flows without 
pressure gradient, oscillating retarded free stream flows Gnd monotonically time-varying 
flows have all three been studied for a wide frequency range. It was found that to the 
lowest frequency considered, the lower frequency hound being determined by economic 
consiaerations (machine time), there weresign if icant vr teady effects or. the turbulent 
boundary layer. 
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The expected phase shifts in wall shear stress and displacement thickness were found, 
as observed by other investigators. Besides the indirect effect of unsteadiness via the 
pressure gradient, the turbulent boundary layer is apparently affected airectly by the 
unsteadiness. This effect is clearly seen in flows approaching zero wall shear stress 
conditions. Unsteady turbulent boundary layers were able to penetrate deeper into 
the adverse pressure gradient whether or not the pressure gradient was alleviated by 
unsteady free stream effects (i.e., onset of zero wall shear stress was delayed in either 
case compared to quasi-steady turbulent boundary layers). For some flow cases, the 
differences were quite large at high frequencies. The unsteady displacement thickness 
and wall shear stress approached the zero wall shear stress condition quite differently 
from the approach of the corresponding quasi-steady values to separation (i.e. no 
singular behavior was apparent in the unsteady case). For the oscillatory flow cases 
considered, the quasi-steady wall shear stress values were a good approximation to the 
corresponding unsteady values for the entire frequency range considered, provided the 
flow was nowhere close to zero wall shear stress conditions. For these same frequency 
ranges, the same cannot be said of Hie displacement thickness values, however. The 
displacement thickness was not represented well by the quasi-steady model even at 
relatively low frequencies. 



INTRODUCTION 


Time-dependence is a conspicious feature of the flow over a helicopter rotor in 
translating motion, and it has lorj been suspected that the effects of time-dependence 
might play an important role in the development of the boundary- layer on such o rotor. 
This suspicion was heightened when it was found (ref. 1) that quasi-steady calculations 
of the boundary-layer development indicated gross separation over nearly one-half of the 
blade disc of a rotor which was known, from experiment, to be unstalled. Although the 
pessimistic outcome of these calculations could not unquestionably be attributed to 
unsteadiness in the boundary- layer, it was clear that a study of the effects of time- 
dependence was urgently needed in order to help isolate the important ingredients of 
this area of rotor aerodynamics. 

The study of the unsteady turbulent boundary- layer is, of course, important in a wider 
context too. Dynamic stall affects, not only helicopter rotors, but also the blades of 
turbines and compressors, and the aerodynamic surfaces of aircraft in maneuvering flight. 
Finally, but not of least significance, the subject is of substantial fundamental interest, 
and a greater understanding of it can be expected to assist indirectly in the understanding 
of a much wider range of boundary- layer flows. 

The subject of unsteady boundary layers is in its infancy, and the subject of unsteady 
turbulent boundary layers was virtually untouched until the last three or four years. 

Since then the differential method of Patel and Nash (ref. 2), and the integral method 
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of McDonold and Shamroth (ref. 3), for calculating two-dimensional flows with both 
spacial and temporal variations, have been published. A few other methods have also 
appeared which treat only temporally-varying flows (refs. 4, 5). 

One objective of this work was to develop further the method of Patel and Nash, ir.d 
extend t to infinite-yowed-cylinder flows as well as two-dimensional flows. A few 
calculations have been done using this new capability. The second objective was to 
perform numerical experiments to explore some of the properties of time-dependent 
turbulent boundary layers, especially those which relate to unsteady separation and 
dynamic stall. Within the scope of this latter objective it was considered important to 
try »o determine the range of validity of quasi-steady methods. Time-dependent 
calculations are more difficult to perform ano are more expensive to perform than 
steady ones, and it is helpful to know the point, in terms of increasing unsteadiness of 
the flow, beyond which the added complexity and expense have to be incurred in order 
to obtain useful results. 
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Description of the method 


Nature ot the Flow Considered 

The method is designed to calculate the time-dependent, incompressible, turbulent 
layer on an infinite yawed cyclinder. Cartesian coordinates are placed on the 
developed surface of the cylinder (whose radius of curvature is assumed to be every- 
where large compart J with the local boundary- layer thickness); x is measured normal 
to the generators, y normal to the surface, and z along the generators. 

The velocities in the x-, y-, and z- directions are expressed in the form U + u, V _r v, 
W + w, respectively, where U, V, W, u, v, w are all functions of x, y and t. They 
are not functions of z because the flow is assumed to be invariant in the direction 
parallel to the generators of the cylinder (ref. 6). The components of velocity U (x/Y/O 
V(x, y, t), W(x, y, t) are defined as ensemble averages, takrn over a 
large number of realizations of the same basic flow, or successive flows with the same 
time-history and the same boundary conditions, e.g. successive cycles of a stable 
oscillatory flow. The components u, v, w represent the random fluctuations about U, 

V, W, and by implication the ensemble averages of u, v, w are identically zero. 

Governing Equations 

The equations of n. tion for a flow of the type described above can be derived from the 
time-dependent Navier- Stokes equations by replacing the three velocity components by 
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U + u, V+v, W+w (as defined above) and forming the ensemble average of tr.e 
resulting set. Introduction of the boundary- layer approx'r. rations, and restriction to 
flows with zero z- derivatives, then yields the two momentum equations for l’ and W 

dij + u £jj_ + v 111 + i_ /_) = o 

at bx by p bx by \ / 

<JW + U dw + V + JL L- ) ~ 0 

at ax ay ay v w / 

together with the statement that the pressure is constant through the boundary layer. 

In the two momentum equations, uv, vw are the ensemble averages of the products 
uv, vw, respectively, and hove essentially the some meaning as the Reynolds stresses 
appearing in the steady turbulent boundary- layer equations. Thus the shear stresses 
T T are given by T and T »-p ™ 

The continuity equation retains its usual farm 

a_u + av_ =0 
a * by 

The above equations are identical to the ones that would have been derived if U, V, 

W had been specified as time averages. However, in order to define time averages, it 
is necessary to make the stipulation that the time scale of the turbulent motion is short 
compared with the time scale of the motion as a whole (ref. 3). A number of flows of 
practical interest involve unsteadiness of the boundary conditions which is too rapid 
for an adequate distinction to be made, particularly when due attention is paid to the 
importance of the large-eddy motions in the boundary layer. These motions can have 
a time scale comparable to the time taken for a particle to be convected a distance 
many times the boundary- layer thickness. 
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To state, as we hove done here, that the equations involve ensemble-overage, rather 
than time-average, quantities, does not mean that it will not be necessary, in due 
course, to make the assumption '\ot ensemble averages can be approximated by time- 
averages. The important point to be made, at this stage, is that the equations are 
mos. secure when viewed as ensemble-average equations. 

The turbulent shear stresses are assumed to behave according to the rate equations used 
in the method for steady three-dimensional flows (refs. 6, 7), but with the convective 
derivative extended to include the time derivative: 
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where $ x '$ z art given by Equation (12) of Reference 7, which leads to co-directionality 
of the shear stress "vector" and the rate -of-s train vector. The assumption of co- 
directionality can be relaxed, if desired, by putting 

$ =$ = 0 . 

x z 

The functions a j, a and L are assumed to be the same functions as apoear in the steady- 
now method (References 6, 7); and L are shown in Figure 1, and Oj = 0.15. This 
assumption implies (a) that the ensemble -average quantities uv, vw can be approximated 
by the time-average Reynolds stresses, and (b) that time-dependence in the mean flow 
has a negligible effect on the structure of the turbulence*. Neither of these hypotheses 
has been verified experimentally, and clearly there is an urgent need to do so. Ft* the 
present, these hypotheses have to be made in order to proceed at all, but it must be 
stressed that some of the conclusions from this study may reed revision if later experi- 
mental work shows that certain underlying assumptions in the method are seriously in 
error. 


Method of Solution 

The governing equations consist of two momentum equations, for U and W, the continuity 
equation, and two rote equations for uv and vw. This set of equations is h}.' >rbolic, and 
is integrated in a three-dimensional domain (x, y, t) by means of an explicit finite- 
difference method.. The calculation scheme is illustrated in Figure 2. The overall 
scheme is similar to that used by Patel and Nash (Reference 2); values of cill the 

*The turbulent energy equation does allow for inertia in transfer of information (energy 
from mean flow to eddy to heat), but the mechanism itself is not a funct : on of the 
unsteadiness. 
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dependent variables are calculated over successive planes, t — constant, and the 
calculation advances in the t-direction . In tr.e present method, however, the dependent 
variables include W and vw in addition to U, V and uv . In addition, the present 

method involves o more sophisticated numerical scheme (Reference 2). and a more 
flexible integration domain; the height of the domain is allowed to vary with x as 
well as t (Figure 3). 

Boundary conditions for the calculation consist of: 

o initial profiles of U, W, uv, vw versus y for t - 0 and all x, 

o initial profiles of U, W, uv, vw versus y for x * 0 and all t, 

o dU = dW = duv = 3vw = 0, and 

dy dy dy d> 

uv - vw = 0 

at y r 1 .258 , 1 .255 being the local height of the 
integration domain, which varies with bath x and t. 
o appropriate boundary conditions at y r 0. 

The wall boundary condition is handled in the some way as in the steady-flow method 
(refs. 6, 7). The outer numerical calculation is matched to an approximate solution, 
based on the low of the wall, at y /$ = 0.05. 

The numerical scheme for the outer layers, 0.05 ^ y/ 5 ^ 1 .25, is equivalent 
to the one-step staggered-mesh scheme described in Reference 8 and is identical to the 
scheme used in the steady-flow method (refs. 6, 7). 
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Five-point central differences are used for the y-derivatives, and either 2- or 3-point 
backward differences for the x-derivatives; backward differences have to be used for 
d/dx, if the U-component of velocity is positive, to avoid the upstream propagation 
of information. 

The maximum permissible t-step is dictated by stability considerations. In most situations 
the t-step is roughly proportional to the smallest value of 5 at the particular t-station 
usually 5 at x = 0. As a result, the computation time varies inversely as the thickness 
of the upstream boundary layer; the smaller the value of 5 at x = 0, the longer the 
run. Figure 4 shows some typical run times on a Univac 1106 computer. For most of 

the calculations, 8 n was taken to be 0.00444, since this value is used in Ref. 2, 

x= 0 

giving a typical computation of c&out 5 minutes. 

For initial boundary- layer thicknesses much smaller than this, the run times become 

excessi.e, and a scheme was developed for segmenting the integration domain so as to 

reduce them. If 5 increases substantially from x = 0 to x = 1 (as it generally does 

if 5^ _ q is small), the computation is inefficient if a constant t-step, determined 

from 5 x q, is used for all x. A continuously variable t-step would be awkward to 

work with, but it is possible to segment the calculation and use a different (constant) 

step-length in each segment. In the scheme tried, segments were chosen such that the 

step-length increased by a factor of 3 from each segment to the next. The factor of 3 

meant that the integration mesh in any one segment merged smoothly into that in the 

adjacent one upstream Figure 4 shows that the scheme was successful in reducing the 

longest computer times, but it could not reduce the shorter ones because the growth of 

the boundary- layer thickness was usually insufficient to permit more than one or two segments. 
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Comparison with Forward Marching Method 


To qualify the accuracy of the computer program, sample calculations were carried out 
for both two-dimeru.jnal and swept infinite cylinder flow fields in the time -re taxation 
mode, in which the boundary conditions are held constant and a steady-state solution is 
approached asymptotically for "large” times. In practice, "large times" turned out lo 
be about 2.5 time units, where a time unit is the time required for a fluid particle 
traveling with some defined characteristic velocity to move a distance of one chord 
length. The results were compared ro solutions obtained by the more conventional 
forward-marching procedure and the correlation is at lea. ' as good as that obtained 
previously by Patel and Nash (Reference 2). 

Figures 5 and 6 show the results of a comparison between the two methods for the case 

of flow over an infinite yawed cyi ; nder, where the chordwise pressure gradient ^p/d x j 

is favorable over the forward part of the cylinder and adverse over the rear. In Figures 

5 and 6, all lengths (x, y, 8*) are made dimensionless by division by the chord-length 

of the cylinder, velocities by division by Q eo (the free -stream velocity at infinity), 

2 

and shear stresses by division by oQq q . Figure 5 shows the predicted chordwise 
variations of 5* T w and T w , The greatest discrepancies between the two solutions 
are 1 .2% in T Wx (at x = 0.26) and 3.0% in 5*(at x = 1.0). Figure 6 shows the 
predicted chordwise and spanwise velocity profiles at x - 1, and it is evident, again, 
that the two methods produce n«jrly identical solutions. 

Based on these results, and similar results for two-dimensional flows, it w as concluded 
that the computer program was performing satisfactorily for flows in which the boundary 
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: conditions are held constant with time. There was no reason to suppose that the same 
would not be true also for flows with time-varying boundary conditions, and it was 
decided to proceed to the calculation of fully unsteady turbulent boundary layers. 


COMPUTATIONAL EXPERIMENTS 

The properties of time-dependent turbulent boundary layers were studied by carrying out 
computational experiments using the calculation method described in the preceding 
sections. Five flow situations were selected, designated Flows A through E, and 
calculations were done for each flow for a range of relevant parameters. Most of the 
calculati . .s were for two-dimensional flows, although an infinite-yawed-cylinder 
version of Flow A was examined and the results are included here. Flows A and B were 
oscillate^ flows; Flows C and D were monotonically time-varying, while in Flow E 
the external flow was first varied monotonically and subsequently held steady allowing 
a relaxation of the boundary layer towards equilibrium conditions. 

The calculations are discussed in detail below, and the results are presented in Figures 

7 through 27. in these figures, cs before, lengths are made dimensionless by division 

by thr length of the plate on which the flow is developing (unity), velocities by division 

by Qe 0 (which is equal to ^e Q in the case of two-dimensional flows), and shear stresses 

2 Z 

by division by pQe 0 (or p Ue 0 in two dimensions). Time is made dimensionless 
by multiplication by U eo and division by the length of the plate. For the oscillatory 
flows a reduced frequency, u>, is defined as 
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o> = 2ir /? 

where P is the period of the motion in time units. If id - 2 tr , a complete cycle takes 
place in the time taken for a fluid particle, moving with velocity U 6o / to be convect- 
ed the length of the plate. In the monotonically time-varying flows, <u is retained as a 
measure of the rate of distortion of the flow, i.e. 1/a > is a characteristic time of the 
motion. 


Oscillating Flows Over Flat Plates 

The flow A is defined as a free stream flow oscillating harmonically, according to the 
equation 

Ue 

= 1 + A sin a> t, 

U e 

e o 

over a flat plate. The amplitude. A, and frequence, w, were varied over a reasonable 
range to determine their effects on the response of the turbulent boundary layer charac- 
teristics to this oscillatory flow field. The initial conditions for this calculation were 
taken to be the asymptotic steady stote solution obtained as time approaches infinity 

(t = 2.5 is sufficient) for A = 0. The boundary conditions for all time at the entrance 

2 

station for the boundary layer (x = 0) were set by making 5 0 - (U /U )”’ and 

e e 0 

scaling the standard initial profiles accordingly. This condition, in fact, is just the 
requirement that the entrance conditions respond instantaneously to the edge velocity 
at any instant of time, i.e. quasi-steady conditions, ’.he Reynolds number based on 
chord length (unity) was taken to be 10^ and the initial boundary layer thickness was 
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.00444. The computer program was allowed to run a few time units (t = 2 seems to be 
sufficient) in order to let transients die out and the output was generated at several 
chord petitions for several periods. Quasi -stead/ solutions were generated by selecting 
tut = constant and utilizing the forward marchiry numerical method for steady-state 
flows. 

Figure 7 shows the results for a ten-fold frequency rrsige at an amplitude. A, of 12.5% 

and Figure 8 shows the same frequency range at an amplitude of 50%. From Figure 7, 

it is seen that the unsteady wall shear stress anticipates the quasi-steady values; i.e. 

there is a phase lead between wall shear stress and the external velocity, U . Similar 

e 

effects have been predicted for oscillatory laminar boundary layers (Reference 9), and 
for perturbed turbulent boundary layers (Reference 3). 

In the laminar case the phase lead increases with frequency, up to a maximum value of 
n/4. In the present results, however, the phase lead is considerably smaller than this 
(see Figure 9). The phase leads predicted here are also smaller than the values calculated 
by McDonald and Shamroth (Reference 3). Their results may have been compromised by 
inadequacies in he velocity-profile model used in their integral method, and the dis- 
agreement may not put the validity of the present results in question. 

At low frequencies, a phase lead is apparent in 8 * also; this phase lead decreases with 
increasing frequency, and at high frequencies the unsteady 8 * actually lags behing the 
quasi-steady values. The average level of 5*, for the unsteady flow, increases w ith 
both increasing frequency and increasing amplitude, until the average unsteady 5* is 
more than twice the average quasi-steady 8* when Aj= 0.5 and uu =15.7 (Figure 8). 
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At the 50% amplitude change (Figure 8) the phase lead in the unsteady wall shear stress 
appears to decrease more with increasing frequency than at the 12.5% amplitude change. 
A peculiar skewing appears in the 5* curve at the higher amplitude although the shear 
stress curves still remain symmetric about their mean values. Apparently, for this case 
the various states the turbulent boundary layer passes through at each instant of time due 
to a change in the external velocity are not reversible when the external velocity change 
is reversed. This situation occurs whenever the condition of zero wall shear stress is 
approached as is the case in Figure 8. The 5* variation for the high frequency case in 
Figure 8 can be approximately represented by a model in which the velocity profile moves 
in lock-step with the external flow. Such a model would be valid as a) -*“00, when 
dU/dt would dominate the other terms in the equations of motion. This "high-frequency 
approximation" allows the statement U 8*= constant. This result is shown in Figure 8, 
the constant having been determined for the best curve fit; it works equally well for the 
lower-amplitude case shown in Figure 7. Whereas the displacement thickness is well 
represented by the high-frequency approximation, the wall shear stress continues to be 
well represented by the quasi-steady approximation, and it is remarkable that it works 
well over such a wide range of both frequency and amplitude. 

Figures 9 through 12 show a comparison of the turbulent results with the larinar results 
of Lighthill (Ref. 9). In Figure 9, a comparison is given of the laminar wall shear stress 
phase lead with the turbulent values versus the reduced frequency, uu. The turbulent wall 
shear stress phase lead is much reduced from the laminar values. Figure 10 shows a com- 
parison of the laminar shear stress with the turbulent values versus reduced frequency, tu. 
Three curves are diown for each case corresponding to three values of the argument tut. 

The values for the high frequency approximation in the laminar case when u)t = rr or 
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3tt/2 are identical and consequently are shown as only a dashed curve. The turbulent 
boundary layer results display the same trends as in the laminar case. The Figures 1 1 and 
12 give graphs of the displacement thickness values analogous to the Figures 9 and 10. 
Once again it is seen that the turbulent boundary layer results display the same trends 
as in the laminar case. 

A variation of the flow A is that of a swept infinite flat plate immersed in a flow oscillating 
harmonically in the chordwise direct^n only, according to the equations 



where Q e is the velocity magnitude at t= 0 and A is the sweep angle. Results of this 
calculation are shown in Figure 13 for the 50% amplitude change. The curves for 8 * and 
(t w ) x have the same characteristics as in ihe two-dimensional case with the exception of 
the skewing of the chordwise wall shear stress at the higher frequency value. This result 
is probably accentuated by the closeness of the chordwise wall shear stress to zero. 
However, the spanwise wall shear stress, in contrast to the chordwise wall shear stress, 
lags the quasi-steady values and this phase lag apparently increases with increasing 
frequency. 
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Oscillating Retarded Flows 


The flow B is an oscillating retarded flow field defined by the equation 

^e = 1 + (A^ ■+• A^ sin a>t)x, 

IT" 

e 

o 

where Aq ■ -0.2. Calculations were made for flow B in the same fashion as for flow A 

with the exception that the initial conditions were taken to be the steady state solution 
obtained as time approaches nfinity (about 2.5) with Aj = 0. Figure 14 shows the 
results for an eight-fold frequency range at the trailing edge when the amplitude change, 
Aj/(1 + Aq), is 12.5% and Figure 15 shows the results for a ten-fold frequency range 
at the trailing edge where the amplitude change is 50%. 

From Figure 14, in contrast to the oscillating flow on a flat plate, flow B gives an unsteady 
wall shear stress and displacement thickness which lag their quasi-steady values. Further- 
more, this phase lag decreases with increasing frequency. It appears that the average 
unsteady 8 * is about the same as the average quasi-stead/ 8* for the entire frequency 
range considered. At an amplitude of 12.5% the boundary layer does not come close 
to zero T w anywhere during the cycle. At the higher amplitude of 50%, however, T w 
becomes very small near the ’kottom" of the cycle. Moreover, and this is most interesting, 
the quasi-steady boundary layer would actually separate during this part of the cycle 
(Figure 15). Thus the unsteady boundary layer manages to remain attached ( T w > 0) 
throughout the cycle, even though the maximum instantaneous retardation would be 
great enough to cause separation if it were maintained for a long time. It is significant 
that, in the present case, the effect of time-dependence is to delay the onset of 
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vanishing t whereas the opposite occurs in the case of the oscillating flat plate (Flow A). 
In that flow, the wall shear stress, at the bottom of the cycle, would eventually decrease 
to zero if the amplitude were increased sufficiently; on the other hand, quasi-steady 
calculations would never indicate separation unless the amplitude were so large as tc 
stagnate the external flow. 

To return to Figure 15, over the part of the cycle where the quasi-steady calculation 
indicates separation, the quasi-steady values of 8* are very large. The unsteady values, 
in contrast, are well behaved throughout the cycle. Figure 15 shows phase lags similar 
to the lower amplitude case in Figure 14 Again, if points close to zero wall shear stress 
conditions are excluded, it is remarkable how good an approximation the quasi-steady 
wall shear stress values are to the unsteady values in a retarded fiow field for a wide 
range of both frequency and amplitude. The results of the high frequency approximation 
for 5* (U e S* = constant) are also shown in Figure 14, and as in Figure 8 the agreement 
with the calculated results is excellent. 


Monotonicaily Time - Varying Flows 


Computational experiments were performed for three flows in which the external velocity 
varied monotonicaily with time, instead of sinusoidally as in Flows A and B. The three 
flows, were two-dimensional, and were defined as follows: 


Flow C: 




a>tx 
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Flow D: 


Flow E: 



U 


e = 1 

U 


e 

o 



+ u>t(1-x) 


+ cut(l-x) for 0 ^ t ^ tp 
+ a>tj(l-x) fo r t > tj 


Only positive values of w and f were considered, and so all three flows involved 
decelerating external flows. In Flow C t! , external velocity distribution "pivoted" 
about the leading edge of the plate (x * 0), and this flow might be regarded as an 
idealized model of the flow in a channel with one wall pivoted and rotating outwards. 


In Flows 0 and E the external velocity distribution pivoted about the trailing edge of 
the plate (x * 1), and these flows might be regarded as idealized models of the flow 
over an airfoil pitching continuously (Flow D), or pitching up to some positive angle of 
incidence and then being held still (Flow E). 

The main objective was to determine the onset of zero wall shear stress and to examine 
the variation of its point of onset with the rate of distortion of the external flow, o > . 
Small values of o> indicate slow distortion of the external flow, and large values indicate 
rapid distortion. As o> — »-0, steady-flow conditions are approached, and fini;^ values 
of tut ^i.e. finite values of dU ^/dxjare reached only for "long" times: t — ^00. 
Steady-f'.ow calculations, done for comparison with the time-dependent ones, could b<» 
performed either by the time-relaxation process or by a forward-marching procedure. 

In most cases the latter method was used for reasons of machine economy; the run-time 
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necessary to reach a steady-state solution in a time-relaxation calculation becomes long 
because of the long convection times associated with the small velocities in the flow 
close to separation. 

The second objective was to examine the general character of the solut * the hope of 
shedding some light on the behavior of the unsteady turbulent boundary approaching 
a point of zero wall shear stress. 

The calculations for all three flows were done for a Reynolds number of 10^ based on U 6q 
and the length of the plate (unity). The initial conditions at t = 0 corresponded to steady, 
flat-plate flow with U e = U e{) and on upstream boundary- layer thickness, jt x — 0 of 
0.00444. For Flow C, the external velocity at x = 0 was constant for all time, and the 
assumption of cjnstant boundary- layer thickness at that point seemed reasonable and this 
was made. However, in Flow D, U e at x = 0 is not constant but increases with time, and 
the problem of constructing reasonable initial conditions there presented a matter of some 
concern. It was decided to try to determine the sensitivity of the solution to changes of 
the initial conditions at x = 0. Calculations were dona for a few representative cases 
with initial conditions as follows: 

(a) flat-plate velocity and shear-stress profiles with constant boundary layer 
thickness (° 0.00444), 

(b) flat-plate velocity and shear-stress profiles with boundary layer thickness 

-0 2 

decreased as ( U /U ) * to account for the change in local Reynolds 

6 e o 

number. 
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(c) velocity profiles, with boundary layer thickness of ' '0444, generated by 

add : ng a constant increment of velocity (- 'J - U ) at all values of y, to 

e e o 

approximate conditions appropriate to tu —00; the hcor stress profiles for 
this case were identical to those in (a) except in the immediate vicinity of 
y ■ 0 (where f w *-••00). 

In all the cases run, vanishing wall shear stress was predicted to >ccur, at x = 1, at 
values of uut which were within 1% of each other for the three sets of starting conditions. 
The calculations indicated larger boundary-layer thickness (up to 36% larger) when (c) 
was used, compared to (a) end (b). However, these larger values of 8 occurred at small 
values of y. : and the differences were much smaller near x = 1 . Nearly identical results 
were obtained with (a) and (b); the predicted values of 8 were the iame to within 1% 
for all values of x greater than 0.02. 

Thus it appeared that the solution was fairly insensitive to the upstream starting conditions, 
and that generality of the conclusions from the present work was unlikely to be compro- 
mised by an unfortunate choice. It was decided that the assumption of constant 8 was 
reasonable, and this was made in all the subsequent calculations. 

Figure 16 shows some results for Flows C and D. The calculations were run until the 

wall shear stress, t w , just reached zero at x = 1, and the value of out at th is nstant was 

3U e 

recorded; uut is equal to the instantaneous velocity gradient . The ratio of 

3 X 

(il)t) T ^ _ o to its value when iu-*-0 is plotted as a function of u) in the left-hand diagram. 
The results show that, for both rlow C and Flow D, t w vanishes, at x ~ 1, ot a higher 
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value of u)t than in quasi-steady flow; in other words, there is a delay in the orv>*?t of 
vanishing t w • In Flow C, (uut) T _ r» increases to c value some 50% greater than that 
for quasi-steady separation, and then remains approximately constant for ui > 0.6. In 
Flow D, the value of (uut) Tw — q increases continuously with increasing m, and there is 
a substantial delay in vonishing T w for large values of uu. For example, v 'th uj = 1C, 
(ujt) Tw = o ' s some five times r’ r eater than it is when at— ^0 (quasi-steady). At tnev? high 
values of utt, the streamwise gradients are so large that the validity of the boundary- 
layer approximations becomes questionable. For this reason, the results for Flow D, 
with at > 8, should be treated with caution, and they are shown by dashed curves in 
Figure 16. The trend suggests that, for high enough values of uc vanishing well shear 
stress might neve r occur. 

One of the effects of time-dependence >s to modify the streamwise pressure gradient, 
dp/Sx. For flow over an infinite yawed cylinder, Euler's equation for flow in the 
external stream takes the form 


dU 

e 

dt 


+ U d u e + i dP =Q 

dx p dx 


The equation shows that 
c-i the sign of dU/ dt. 

R = -( 


dp/dx can be either increased or decreased, depending 
The qur.ntity R, where 


dU / 
e 


L i_L\ 

p ' 


is a measure of the contribution of . >e total pressure gradient. 
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Posit! *e values o* R correspond to coses where dp/dx is increased, i .e. becomes 
more adve».«*, negative values to coses where the adve r se pressure gradient is alleviated. 

The r’jht-hand diagram in Figure !6 shows the results for Flows C, D ploted versus R. 

Act* Jly, R »ss function of x, and the valy^s plotted correspond to x = 0.5. Flow 

D generates negative values of R and the delay in the onset of zero t takes place in 

an environment where the adverse pressure gradients are alleviated by the effects of 

time-dependence. On the other hard. Flow C generates positive values of R and, 

while the effects are less marked than in Flow D, there is again a delay in vanishing 

T , and it occurs in the face of an increasingly adverse pressure gradient, ’ndeed, 
w 

at Hie highest value of to considered the value of bp/ bx is some 300 times greater 
than was necessary to cause separation of x = 1 in steady flow . 

The inference is that the augmentation or alleviation of bp / bx G°e$ play a part 
in *he effect of time -dependence on vanishing wall shear stress, but that, at least 
for flows similar to the present ones, it is not the only factor to be considered. The 
fact that a delay in vanishing occurred even when bp I d x was increased 

suggests that time-dependence influences the flow directly as well as via bp / d* • 

It might also be worth mentioning that, regardless of the sign of b U / b t and its 
effect on bp/ b* / the point of vanishing f can move upstream only if b\J/ b t 
in the vicinity of the wall is locally negative. 

In Figures 17, 18/ T and 5*, from the calculations for Flow C, are plotted versus 
x for a low (0.25) and a high (31 .4) value of cu . The results correspond to the instant 
of time, and the instantaneous value of cut, where T just vanishes at x = 1 . The 
results of two quasi-steody calculations are also plotted for comparison. In one, the 
value of tut is the same; in the other, the value of cut is reduced such that separation 
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occurs at x = 1 . With ut the same, separation of the steady boundary layer occurs 
ahead cf x = 1 (vanishing wall shear stress is delayed in the unsteady case)* In order 
to obtain separation at x = 1 , it is necessary to reduce <t,t, i.e. to reduce the stream- 
wise adverse pressure gradient . 

Figures 19, 20 show the corresponding results for Flow D, again for a low and a high 

value of cj . It will be noted, in Figure 20 especially, that the values of for the 

two quasi-steody runs are substantially different from one another particularly ot low 

values of x. This difference is a consequence of the different levels of between 

the two quasi-steady flows; in Figure 20, for example, the two values of cut used in 

the quasi-steady calculations are different by a factor of 9. Plotting t /( pU ), 

w e 

2 

rather than t /( oU ) would have reconciled the results for small x, but would 
w K e 

o 

have obscured important trends near separation. 

All four sets of results (Figure 17 through 20) exhibit essentially Hie same characteristics. 

The steady-flow separation is associated with increasing values of -( ^ x / dx), 

w 

as the separation point is approach, reminiscent of the square-root singularity observed 
in the laminar case. Correspondingly, 8* increases rapidly, again, in a quasi-singular 
manner. In contrast, in the time-dependent cases, the approach to zero wall shear stress 
is more gradual: the value of -( dt / dx) decre~ e e$ rather than increases with x. 

W — 1 

The results suggest that the delay of vanishing is directly linked to the absence of 
the square-root singularity (or something close to it) in the t (x) variation. The 
curves of 8* also depict important differences. Particularly at the higher values of 
to , 8* is smaller and has a markedly smaller slope as goes to zero. In short, 

there is no evidence of singular behavior, and thus the results support the views of 
Sears (Reference 10). It was not possible to continue the calculations beyond the 
point of instantaneously zero wall shear stress, and so the question of the existence 
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of a singularity downstream of this point could not be addressed. Neither could .aU 
culations easily be done for the situation where the point of zero is moving down- 
stream. For this case, the Seers model (Reference 10) of unsteady separation wouid 
predict that the singularity lies ahead of the point of vanishing t . 

The results presented so far, for Flows C. D, all correspond to the instant of time ci 
which = 0 at x = 1 . Further studies were conducted to examine the upstream 
movement of the point of vanishing x w to smaller values of x. 

Figure 21 shows curves of i versus x for increasing time; these are for Flow D with 

w 

oj = 1 .57. As stated before, increases rapidly with time, for low values of x, 

because of the increase of U . When i fell to zero at x = 1, the row of node 

e w 

points of the integration mesh, corresponding to x = 1 , was dropped, and the calcula- 
tion was continued to largervalues of t. 

Additional rows of node points were dropped as the point of zero x^ moved upstream. 
Figure 22 shows the corresponding curves of 5* versus x. It will be noted that, up to 
a time t = 1 .9, £* is independent of time for x less than about 0.2, This result pro- 
vides indirect confirmation that the assumption of constant 5 (and therefore nearly 
constant 5*) at x = 0 was a reasonable one. It will be noted, also, that the curves 
become progressively steeper as the point of zero moves forward. It is not clear, 
at this stage, whether a singularity is being re-established, or whether this effect is 
due to numerical error as the result of there being fewer node points left in the inte- 
gration net. 

Figure 23 illustrates the velocity and shear stress profiles calculated for Flow D (with 
c o = 1 .57) at x = 0.67. The wait shear stress vanishes at x = 0.67 when t = 1 .93, 
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i.e. when cut is approximately 3. Quasi-steady separation would occur ot x = 0.67 

only if i>t were reduced to 1 .55. The profiles corresponding to this latter situation 

ore shewn for comparison. It is interesting to note that the velocity profiles ore similar 

in thickness and shape, the only r» nificant difference being the edge velocity, U ; 

thus plots of versus y would look almost identical . Fu.thermore, ' ihear 

stress profiles ore similar in shape and would fall close io one another if plotted in 

2 

terms of t/( p U ) . 

This similarity in shape is surprising in view of the difference in behavior of <r (x) 

w 

and $*(x) as the point where <■ = 0 is approached. 

The same time-dependent flow (Flow D, with o> - 1 .57) reaches cut = 1 .55 when t is 
approximately one. The profiles corresponding to this condition are shown in Figure 23 
for comparison both with the quasi-steady results, at the some cut, and with the results 
for Flow D at the later time, t = 1 .93. The comparison shows that the two velocity 
profiles with the some edge velocity, U^, are substantially different in thickness and 
shape. Moreover, the maximun. shear stress is quite different. 

Thus it would appear that the similarity between time-dependent and quasi-steady flows 
is closest when they ore compared at the points where *r = 0 in both, even though 
the edge velocities are then different and even though the upstream histories of the two 
flows are radically different. Conversely, development in rhe same external velocity 
distribution appears to be largely irrelevant; it does not ensure >he same upstream history 
and does not ensure even approximate similarity between the fwo flows. 
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Figure 24 shows the variation with time of < and 8* at x = i f tor Flow D with 

w 

Co = 1 .57; <r goes to zero when t - 1 .37. Flow E, the results of which are shown 
w 

for comparison, is ider tical to Flow D for 0 - t ^ 1 , but for t >1 the external flow 

is frozen with cut = 1 .5/. In Flow E, <r goes to zero at t = 1 .47. 

w 

Figure 25 shows the forward movement of the point of zero for Flows D, E with 

(i) = 1.57, and also with a higher value, cu = 31 .4. In Flow D the ext 3 rnal velocity 

gradient increases continuously, and the point of zero moves forward continuously 

at least as far as x = 0.5. The higher the value of cj , the more rapidly this point 

moves forwar \ In the case of Flow E, the point of vanishing moves forward to 

the position of the quasi-steady separation point for the particular value of cut; 

here, around x = 0.67. Thus for t >1 when a>= 1 .57, and for t > 0.05 when cj = 31 .4, 

the external velocity distribution, U (x), is frozen and the movement of the point of 

e 

zero -f is just part of the relaxation towards steady-state conditions. It is inter- 
esting to note that the time scale describing this relaxation process appears not to 
be very dependent on the value of to corresponding to the previous unsteady external 
flow. It is also interesting to note the comparatively long time necessary for the flow 
to settle down to the new steady-state situation, nearly 4 time units, which (at this 
value of cut) is roughly five times the time taken for a fluid particle in the external 
flow to be convected from x = 0 to x = 1 . 

During the relaxation phase of Flow E, while the external velocity distribution is 
frozen, the characteristics of this flow have to change from those appropriate to 
time-varying conditions to those appropriate to steady-state conditions. This 
process is illustrated by Figure 26 which shows the variation of T and 8* with 
x and t for Flow E with to = 1 .57. The shape of the ^(x) curves changes from 
one in which -( 3 x) decreases with increasing x to one where it increases 

with x. At the same time, the curves of 8* versus x become steeper ns the point 
of zero r is approached. 
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Figure 27 shows the corresponding velocity and shear stress profiles. It is interesting to 

note how the velocity profile approaches that for steady-state separation by a combination 

of retardation and thickening which leaves the average velocity gradient, d*j/Sy, nearly 

constant. The maximum shear stress increases far more than the wall shear stress de- 

2 

creases. It increases roughly in proportion to 8 , rather than to 8 which would be 
predicted from mixing-length theory for a boundary layer with a constant dU/dy. This 
departure from a mixing-length type correlation between the shear stress and the velocity 
gradient is to be expected in this flow because of the large streamwise gradients. 
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CONCLUSIONS AND RECOMMENDATIONS 


The method developed here for calculating time-dependent turbulent boundary layers in 
incompressible flow past infinite yawed cylinders is a logical development from the work 
of Nash end Patel on steady three-dimensional turbulent boundary layers. One of the 
space dii"t.nsions of the integration domain is replaced by time, while the three- 
dimensional, ty of the velocity and shear-stress vectors is retained; U, V, W, <r , 

X 2 

are determined in the domain (x, y, t). 

The empirical content of the method, the dissipation length, diffusion function and the 
ratio of the shear stress to the turbulent kinetic energy, are carried over from the steady- 
flow methods. The assumption is made that the original forms of these functions still 
hold, and the additional assumption is made that the ensemble averages in the time- 
dependent equations can be approximated by time-averages. These assumptions have 
not been verified experimentally, and until they are, the conclusions of the present 
study cannot be accepted without reservation. However, the aim has been to establish 
trends rather than definitive values, and it is not unreasonable to expect that these 
trends will survive possible changes in the empirical input which are called for as the 
result of further experimentation. 

It is quite evident, from the results shown, that there are significant unsteady effects on 
the turbulent boundary layer particularly in the regime of approaching zero wall shear 
stress. For the oscillatory cases examined, even for the lowest frequencies calculated, 
there were significant effects on the displacement thickness. A small degree of phase 
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shift was observed between the wall shear stress and the external velocity. However, 
the level of wall shear stress differed little from quasi-steady values, corresponding to 
very low frequency, until zero wall shear stress was approached. 

The effect, on the boundary layer, of time-dependence in the external flow is complex, 
and cannot be represented merely by a change of the streamwise pressure gradient. For 
some flows the onset of vanishing wall shear stress is delayed, compared with the 
equivalent steady flow, not only when the adverse pressure gradient is alleviated by 
the unsteadiness, but even when it is increased. The delay was clearly demonstrated 
by the monotonically time-varying flow examples where the approach to zero wall 
shear stress conditions was examined. Unsteady turbulent boundary layers were able to 
penetrate deeper into adverse pressure gradients (i.e. delayed onset of zero wall shear 
stress) than could be accounted for by quasi-steady calculations. For some flow cases, 
the differences were quite large at high frequencies. Ir. addition, the unsteady dis- 
placement thickness and wall shear stress approached the zero wall shear stress 
condition quite differently from the approach of the corresponding quasi-steady values 
to separation. The results support the contention of Sears (Reference 10) and others that 
the point of vanishing wall shear stress, in an unsteady boundary layer, is not a singular 
point. Furthermore, the results suggest that the fact that it is not a singular point has 
a lot to do with the delay in vanishing wall shear stress mentioned earlier. 

This study did not address the question of separation, end we have been careful to avoid 
any suggestion that the point of vanishing wall shear stress is a separation point in 
unsteady flow. Indeed, if separation, in the sense of flow detachment, is associated 
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with singular behavior of the displacement thickness, the present results indicate that 
separation was always downstream of the point of vanishing wall shear stress for the 
flows examined. Thus it would appear that separation was delayed even more by the 
effects of time-dependence than was the onset of vanishing wall shear stress. There is 
an urgent need for a study of the events surrounding actual separation, i.e. detachment, 
o^ the unsteady turbulent boundary layer. With certain modifications the present 
calculation method could be used to examine the development of the flow between the 
onset of zero wall shear stress and the onset of separation. 

This study has drawn attention fo the strengths and also to the weaknesses of the present 
calculation scheme. The uncertainties associated with the f'ow model have already 
been mentioned, but a few comments are in order concerning the numerical scheme for 
integrating the governing equations. The explicit scheme used here is well tried, 
having been employed in a range of three-dimensional steady boundary- layer calculations 
as well as in the present unsteady-flow method. There are deficiencies in the scheme 
from the standpoint of precision, particularly in the inner-layer calculation; however, 
the numerical precision is probably at least as good as the accuracy of the empirical 
flow model. The major disadvantage of the present scheme lies in the long computarion 
times required. It has not been possible to explore important low frequency ranges 
because the calculations became too expensi- e, nor was it possible to treat flows where 
the boundary- layer thickness at the upstream end of the domain was small. There is an 
urgent need to improve the economy of the present method, and it would seem that this 
can be done only by replacing the explicit numerical scheme by an implicit one which 
makes no demands on step size for stability. The conceptual framework of a suitable 


31 



implicit scheme has already been established. Be that as it may, the existing method 
represents a very useful tool, as it stands, and a number of further investigations c n be 
made even before the implicit method is in an operational state. The tools developed in 
this analysis should be used to investigate the differences between the turbulent boundary 
layer produced by an unsteady free stream but stationary wall and a steady free stream 
but non-stationary wall. Of considerable interest would be the utilization of these 
tools to calculate unsteady turbulent boundary layers which correspond to the experimental 
investigations being conducted by the U.S. Army Air Mobility R&D Laboratory and the 
Naval Postgraduate School in Monterey. The comparison of these experimental data 
with the theoretical values would not only serve to elucidate unsteady turbulent boundary 
layers but would also point out ways in which the mathematical model might be improved. 

Aside from these fundamental investigations, the tools developed in this research program 
are ready to and can be utilized to great advantage in studying the practical problem of 
pitching airfoils and dynamic stall. In particular, an investigation can be carried out 
on the question of lift coefficient overshoot with regard to when and how does separation 
occur. Related problems which can be studied are the effect of frequency on displace- 
ment effect and the resultant effects on pressure distribution. It is entirely possible that 
current tools for calculating the unsteady turbulent boundary layer can shed some light 
on the bursting of the separation bubble. The capability now available for unsteady 
turbulent boundary layers coupled with advanced methods for treating the potential flow 
could well answer the question of how the boundary layer rolls up after separating from 
the pitching airfoil . 
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The analysis and computer program can and should be extended to compressible flow. 

When that is accomplished, the required technology for examining the very complex 

problem of shock-unsteady turbulent boundary layer interaction would be in hand. 

SUMMARY OF CONCLUSIONS 

1 . There are significant unsteady effects on the turbulent boundary layer particularly 
in the regime of approaching zero wall shear stress. 

2. For the retarded free-stream flow case, the effect of time-dependence is to delay 
the onset of vanishing t w whereas the opposite occurs in the case of the oscillating 
flat plate. 

3. Unsteady turbulent boundary layers were able to penetrate deeper into adverse 
pressure gradients whether or not the pressure gradient was alleviated by unsteady 
free stream effects (i.e., onset of zero wall shear stress was delayed in either case 
compared to quasi-steady turbulent boundary layers). 

4. The unsteady displacement thickness and wall shear stress approached the zero 
wall shear stress condition quite differently from the approach of the corresponding 
quasi-steady values to separation. 

5. In contrast to the steady-flow separation, associated with increasing values of 
“(bT w /?x), the time-dependent approach to zero well shear stress is characterized 
by -(dT w /dx) decreasing. 
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6. (f separation/ in the sense of flow detachment, is associated with singulai behavior 
of the displacement thickness, the present results indicate that separation was 
always downstream of the point of vanishing wall shear stress for the flows 
examined. 

7. Regardless of the sign of dl)e/dt and its effect on dp/dx, the point of vanishing 

t w can move upstream only if dU/dt in the vicinity of the wall is locally negative. 

8. For FSw D, the solution was fairly insensitive to the upstream starting conditions 
allowing the assumption of i constant initial boundary layer thickness to be made. 
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Figure 3. Section x - y of the Integration Domain for a Typical Calculation 






Figure 5. Comparison of Time-Relaxation Results with Forward 
Marching Results for a 30 Swept Infinite Cylinder 
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Figure 7. Oscilla 
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Figure 8. Oscillatory 
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Figure 9. Wall Shear Stress Phase Lead Vs. Reduced 
Frequency for Laminar and Turbulent Flow 
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Figure 10. Wall Shear Stress Vs. Reduced 
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Figure 12. Displacement Thickness Vs. Reduced Frequency 
for Laminar and Turbulent Flow 





Quasi "Steady 



Figure 13. Oscillatory Flow Over 30° Swept Flat Plate At High Amplitude 
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Figure 14. Oscillatory Retarded Flow at Low Amplitude 
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® Unsteady, (X) _ = 1 . , cu = .251 , (cut) = 1.218 
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® Steady State, same cut, different (x) 



Figure 19. Wall Shear Stress and Displacement Thickness 
Distribution for Flow D at Low Frequency 
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Figure 20. Wall Shear Stress and Displacement Thickness 
Distribation for Flow D at High Frequency 
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Figure 25. Upstream Movement of the Point of Zero Wall Shear Stress; Flows D and E 
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